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Abstract 

We analyse nonperturbatively signal transmission patterns in Green's functions of 
interacting quantum fields. Quantum field theory is re- formulated in terms of the 
nonlinear quantum-statistical response of the field. This formulation applies equally 
to interacting relativistic fields and nonrelativistic models. Of crucial importance is 
that all causality properties to be expected of a response formulation indeed hold. 
Being by construction equivalent to Schwinger's closed-time-loop formalism, this 
formulation is also shown to be related naturally to both Kubo's linear response and 
Glauber's macroscopic photodetection theories, being a unification of the two with 
generalisation to the nonlinear quantum-statistical response problem. In this paper 
we introduce response formulation of bosons; response reformulation of fermions 
will be subject of a separate paper. 
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1 Introduction 



This paper continues the investigation of response properties of quantum sys- 
tems started in Ref. In [ij we introduced a response formulation of the 
harmonic oscillator and extended it to noninteracting bosonic fields. Here we 
show that response formulation may be further extended to arbitrary inter- 
acting bosonic fields. Response reformulation of fermions will be subject of a 
separate paper. 
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For noninter acting bosons response formulation means description of the 
quantum system in terms of quantum averages of the normally- ordered prod- 
ucts 0, 0, [1, 0] of field operators defined in the presence of external sources. 
In we proved that this description is equivalent to the standard quantum- 
field-theoretical description of the same system within Schwinger's renown 
closed-time- loop formalism [6|. The response formulation and the Schwinger 
formalism are coupled by a one-to-one response substitution in the correspond- 
ing characteristic functionals. 

It would seem that the obvious way of generalising this result to interacting 
fields is replacing the normal ordering of free-field operators by the time- 
normal ordering of Heisenberg operators as introduced by Glauber and Kelly 
and Kleiner [3, [5|. However, we quickly discover that this leads to loss of 
the key property of the response formulation: its equivalence to the standard 
Green-function approach. Without an amendment to the concept of time- 
normal ordering, response substitution for interacting bosons does not exist. 
This makes the whole exercise pointless: recall that our ultimate goal is ex- 
tending the phase-space techniques to relativistic problems. We are therefore 
forced to choose a different approach. We take the response substitution found 
in [l| for noninteracting bosons, postulate it for interacting bosons, and con- 
sider the consequences. This implies introducing a new definition replacing the 
familiar time-normal operator ordering. However, we also show that within the 
optical paradigm (technically, within the approximation of slowly varying am- 
plitudes) this definition coincides with the definition by Glauber and Kelly 
and Kleiner. Interpretation of any of the quantum-optical experiments needs 
not be reconsidered. 

Except being a preparatory work for the phase-space approach to relativistic 
quantum fields, some results of this paper appear to have significance of their 
own. First and foremost, we demonstrate that all quantum properties of in- 
teracting systems may be interpreted in terms of response and self-radiation. 
As was explained in paper one interpretation of our results is proving the 
equivalence between Schwinger's closed-time-loop formalism j6| on the one 
hand, and a certain generalisation of Kubo's and Glauber's approaches com- 
bined on the other. For more details we refer the reader to the introduction of 
Ref. All points made there apply not just to the harmonic oscillator, but 
also to any interacting bosonic quantum system. 

Perhaps the most interesting result of this paper is a fundamental link between 
response and noncomutivity of operators. Indication of this connection naay 
be seen already in Kubo's famous formula for the linear response function [7|, 




where the latter is expressed by the average of the two-time commutator. This 
feature is shown to hold for the full nonlinear quantum-statistical response 
of interacting systems. The assumption that operators commute cancels the 
dependence of the system properties on external sources turning the system 
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into a "pre- assigned quantum source." 



Another interesting result, which plays only a technical role in our analyses but 
seems to be important on its own, is that all response properties of a system 
are contained in the field operator defined without the external source. In 
other words, the information contained in the field operator (and in the intial 
condition) suffices to describe any scattering experiment performed on the 
system, a fact which has never been fully appreciated and which has profound 
consequences for the quantum measurement theory. Again, a known example 
of this is Kubo's formula where the linear response to a source is specified in 
terms of the field operators defined without the source. 

One technical point merits special mentioning here. All usual problems of the 
quantum field theory like adiabatically switching the interaction on and off, 
renormalisations and the like are securely locked away in the assumption that 
Heisenberg field operators and Green's functions "may be defined." In partic- 
ular, we do not make any effort at specifying Green's functions at coinciding 
time arguments because in all practical calculations such specifications emerge 
as a side-effect of a regularisation procedure. This is the case in the relativis- 
tic quantum field theory [sl as well as for simple nonrelativistic models j^. 
Note that, in the latter case, regular isat ions may be applied directly to the 
phase-space equations in the form of specifying a stochastic calculus j9|. A 
similar regularisation procedure will be introduced in our forthcoming papers 
for relativistic phase-space models. 

The paper is organised as follows. In section [2], we introduce the necessary 
quantum-field-theoretical concepts, such as time orderings of operators. Green's 
functions and their characteristic functionals. In section [3|, we reiterate the key 
results of paper [i]. Using these results as leading considerations, we then pro- 
ceed to defining the response formulation of interacting bosons. In section H] 
we prove that the two ways of describing the system, in terms of the field op- 
erator defined with and without the external source, are formally equivalent. 
In section |5l we define the concepts inherent in the response formulation, such 
as time-normal operator ordering, quantum-statistical response functions and 
related characteristic functionals. In particular, we show why the Glauber- 
Kelly-Kleiner definition of the time-normal ordering must be amended, and 
why our amendmend is inessential for quantum optics. In section [6l we de- 
rive formulae relating Green's functions to response functions and vice versa. 
Among other things, we show that our expression for the linear response func- 
tion coincides with Kubo's formula, and give explicit demonstration of the 
link between operator noncommutativity and response. In section [71 we prove 
explicit causality in the response formulation. In appendix \^ we summarise 
the necessary facts regarding separating the frequency-positive and negative 
parts of a function. Finally, in appendix [HI we summarize the results from the 
response formulation of charged bosons, while in the main body of the paper 
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Fig. 1. Schematic diagram of a response experiment. The source and the detector 
are coupled to a quantum field (dashed lines with arrows), e.g., the electromag- 
netic field. This field also interacts with some matter and/or environment which 
shapes its observable properties. The dotted rectangle encircles the part of the sys- 
tem (field+matter+source) we treat explicitly while a detector is implicit in our 
considerations. 

we only treat neutral bosons. 



2 Statement of the problem 

2.1 Neutral bosons 

2.1.1 Quantum dynamics in the presence of a source 

The Gedanken experiment we have in mind is illustrated in Fig. [H It involves 
two classical devices, the source and the detector, coupled to a quantum field 
described by a Hermitian 4- vector field operator Q,^{r, t). Formally, we analyse 
a quantum- statistical response of the quantum field (cf. Fig. [T]) by considering 
a system with the generic Hamiltonian 

n{t)=no{t) + Y. I d'rfir,t)Q,{r,t)=no{t)+3{t)Q{t). (i) 

Here, the external source j{t) is a given real c- number current, and we have 
introduced a general shorthand notation, 

<t)Yit) =Y.j d'rx'^ir, t)Y,ir, t) = E / d'rx,{r, t)Y''{r, t), (2) 
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where the integration extends the whole 3-dimensional space, cf. endnote 
Shorthand ([2]) apphes by default throughout the paper. 
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By definition, Q^{r,t) and 'Ho{t) are the field operator and the Hamiltonian 
in the Heisenberg picture with j = 0: 



(3) 



The field operator and the full Hamiltonian in the Heisenberg picture with 
j 7^ will be denoted as Qj^{r,t) and Hjit): 



(4) 



To assign formal meaning to the quanties used in Eqs. ([3]) and (jlj), we assume 
the existence of the Schrodinger-picture Hamiltonian, 



ns{t) = nos{t)+j{t)Q, 



(5) 



where Tiosit) includes the free- field Hamiltonian as well as all Hamiltonians 
describing the matter and environment in Fig. [1] with corresponding interac- 
tions. The interaction of the field with the external source is described by the 
term j(t)Qs, 



(6) 



where Qsfi{r) is the field operator in the Schrodinger picture (assumed to be 
time-independent). Then 



Q^,{r,t)=lJ\t)Qs,{r)lJ{t), 
Mt)=Ul{t)H,sit)Uo{t), 



(7) 



where the evolution operators U{t) and Uoit) related to iisif) and Ti.os{t), 
respectively, obey the Schrodinger equations. 

It is assumed that U{t) 1 ast ^ — oo; this applies to all evolution operators 
defined here as well as in section |H 

The fact that all nonlinear interactions enter through Tiosit) in no way pre- 
vents us from treating this term as the unperturbed Hamiltonian and the linear 
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term j(t)Qs as interaction. Following the standard time dependent perturba- 
tion techniques, we define the interaction-picture evolution operator Uj{t) as 
the "difference" between W(t) and ^0(^)5 



(9) 



The evolution operator Uj{t) obeys the equation (using shorthand 



(10) 



solved by the T-exponent, 



Ujit) = T+ exp 



hJ- 



dt'3{t')Q{t') 



'111 



We can then relate Q.j^{r,t) to Q^{r,t) directly, 



Q,-^(r,t)=Wj(t)Q^(r,t)W,(t). 



'12) 



For purposes of this paper, it suffices to assume that the Heisenberg opera- 
tors Q^(r,t) "in the absence of the source" and Qj^(r,t) "in the presence 
of the source" are known. Equations (ITT]) and 0121) show that Qj^{r, t) may 
be expressed by Q^{r,t)^ so that the only assumption we really need is that 
the operator Qfj,{r,t) is known. All other details, including the actual phys- 
ical nature of the operator Q^{r,t), are irrelevant. Our choice of a 4-vector 
field suggests quantum electrodynamics, however, we do not impose gauge 
conditions nor symmetry nor any transformation properties. What we do here 
equally applies to a full relativistic bosonic field and a single resonator mode. 
The index fi can in fact be anything, and run over any arbitrary set of indices, 
modes or polarisations. 

Furthermore, we do not make any assumptions regarding the Hamiltonian 
7ios(^) except its existence. It may include arbitrary interactions with other 
quantum fields (cf. "Matter+Environment" in Fig. [1]). These interactions may 
be explicit, or accounted for phenomenologically as nonlinearities, or treated as 
heat-baths, or as any combination of these. Nothing precludes them from being 
explicitly time-dependent. In other words, our analyses apply to arbitrary open 
nonlinear systems. We do not assume that closed evolution equations of any 
kind exist for the field operators. This frees us of all physical restrictions, be it 
energy conservation or Markovian nature of the evolution. The only applicable 
condition is conservation of probability. 
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2.1.2 Closed-time-loop formalism 



The minimal approach capable of accommodating this kind of generality is 
Schwinger's famous closed-time-loop formalism 0]. Our quantities of interest 



are the Schwinger-Perel-Keldysh-style [10| double-time-ordered averages of the 
field operators 



T-Q,AriM) ■ ■ ■ Q,Jr^,UT+Q,'y,,t[) ■ ■ ■ Q.^S^'uX) )■ (13) 

Time-ordering (T_) puts operators in order of increasing (decreasing) time 
arguments. By definition, bosonic operators under the ordering sign commute, 
and, if ti < t2 < • • ■ < tm, 

T-Q{h) ■ ■ ■ Q{tm) = Q{tl) ■ ■ ■ Q{tm-l)Q{tm), 

T+Qih) ■ ■ ■ Qit^) = Q{t^)Q{t^.i) ■ ■ ■ Q(ti), 

where all other arguments of the field operator are omitted for brevity. Equiv- 
alence of the double-time ordering with the closed-time-loop ordering was 
discussed in section 2.2 of Ref. fl|. The quantum averaging in (IT^ is defined 
in the standard way with respect to the Heisenberg density matrix 

(•■■) = Trpo(---)- (15) 

This definitions implies the Heisenberg picture with respect to Hamiltonian 
7^0- For reasons explained in the introduction, we do not make any effort at 
specifying Green's functions at coinciding time arguments. 

The whole assemblage of Green's functions fll3p is conveniently handled through 
their generating, or characteristic, functional r^^): 



(-i)"^i" 5™+"$(r7_,r7+) 



, (16) 

r)-=r]+=0 



$(?7_,?7+) = (T.exp 



1+ E 



m+n>l 



mini 



dtri4t)Q{t) 

m+n 



T+ exp 



-I / dt7]+{t)Q{t) 



dt 



v-{ti) ■ ■ ■ v-{tm)v+{t'i) ■ ■■v+iQ 



X (r_Q(tO ■ ■ • Qit^)T^Q{t[) . . . Q{Q) . (17) 



The functional arguments rjt^ir, t) here are arbitrary smooth 4- vector c-number 
functions, and the symbol (/ dt)"^~^^ denotes integration over all time argu- 
ments. Omitting integration limits means that the time integrations are from 
minus to plus infinity, cf. also endnote [13|. Equation (JTTll implies shorthand 
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([2]); it also gives an example of a "loose" usage of this shorthand, with the 
corresponding pairs (e.g., ri^{ti) and Q(ti)) separated by other factors. Defi- 
nitions identical to (fT5]) - (IT7|) may also be written with operator Q.j^{r,t) in 
place of Q^(r,t). This introduces another functional 



exp 



dtri-{t)Qj{t) T+exp - i J dtr]+{t)Qj{t) 



:i8) 



so that $(?]_, = $(r7_, r7+; 0). For simplicity, we assume that 

Jit) = 0, \t\ > Wx, (19) 

where tmax is some large parameter, and that the Heisenberg pictures with 
and without the external source coincide for t < — tmax- This allows us to use 
the same averaging in (ITTl) and (ITS!) . 



Shorthand ([2D plays in our analyses much more important a role than merely 
saving space. It effectively reduces the problem of the quantum field to the 
generic case of nonlinear quantum oscillator interacting with environment, 
with Q{t) being the oscillator coordinate. In turn, the analyses introduced 
with the nonlinear oscillator suffice because our formal manipulations will 
always apply to the time argument of the field operator leaving the rest of the 
field arguments (labels) alone. All formulae in the analyses below belong to one 
of two groups. Firstly, relations for characteristic functionals like (ITTl) apply 
in the field case directly by virtue of definition ([2]). Secondly, relations for 
operators and quantum averages like (|T^ are generalised by simply restoring 
all labels 

t^^,r,t, dt d^rdt. (20) 

•' /I 

We shall consistently refer to Q{t) as "field operator" so as to emphasise 
that our results apply not just to the nonlinear oscillator but to an arbitrary 
bosonic field. However, treating time as variable and space as label makes our 
viewpoint essentially nonrelativistic. All formulae in this paper apply, strictly 
speaking, in a given reference frame. 



2.2 Charged fields 



A charged field is described by introducing a Hermitian-conjugated pair J-'^{r, t), jFt(r, t) 
of Heisenberg operators where 



:Fl{r,t) 



:F^{r,t) 



(21) 
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Above definitions for neutral bosons can be generalised to a charged bosonic 
field by simply redefining the linear forms of field operators: 



j(t)Q(t) = E / rfV[/*(r,t)J-^(r,t)+ J-t(r,t)/(r,t) 



(22) 
(23) 



cf. endnote 3]. This equals treating as a two-component field, or, 

putting it another way, introducing an additional "label of complexity" dis- 
tinguishing T from P . Replacing Q by a pair of operators leads to doubling 
all functional arguments. Functional f[T7l) is replaced by 



T_ exp <i dt 



X T+ exp \—i dt 



■ (24) 



The definition of Green's functions f[T^ is amended accordingly. Given the 
triviality of generalisation to charged bosons, and to prevent our formulae from 
becoming even bulkier, in the main body of the paper we only consider neutral 
bosons. Summary of the formulae for charged bosons is given in appendix [Bl 

The case of fermions (which are always charged) may be reduced to the case of 
charged bosons by redefining operator forms fl2^ and (1231) with the r/s and js 
being g-numhers in place of c-numbers. G- numbers, or Grassmann variables, 
are often called classical anticommuting quantities, meaning that g-numbers 
mutually anticommute as c-numbers mutually commute. However, introducing 
Grassmann currents into Hamiltonians requires a number of modifications of 
the conventional Hilbert space techniques. Another subtlety is that definitions 
of Green's functions for fermionic fields involve sign conventions. Response 
properties of fermions will therefore be a subject of a separate paper. 



3 Harmonic oscillator revisited 



Here we briefly summarise the response formulation of the harmonic oscillator 
introduced in Ref. For the oscillator, Q{t) q(t) and Qjit) —>■ qjit), 
where q{t) and qj{t) are the Heisenberg-picture displacement operators defined 
with and without the source. Definitions f[T7|) and f[T5]) then apply with the 
replacements Q.{t) q(t) and Qj{t) qjit), and do not assume shorthand 
([2]). We do not introduce any special notation for the functionals $(?]_, r^^) 
and $(r7_, j) for the harmonic oscillator. Instead, the fact that a particular 
relation applies only in the oscillator case will always be stated explicitly. In 
this section, all formulae apply to the harmonic oscillator. 
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The response formulation of the oscillator is based on the functional $r(?7; j) 
defined as 



$R(r/;j) = ^:exp J dtr]{t)q,{ty}j . (25) 

Operators q{t) and qj{t) differ by a c-number displacement under the influence 
of the external source, 

q,{t) = q{t) + q,{t), q,{t) = J dt' D^{t - t')j{t'), (26) 

where -DR(t) is, up to the overall sign, the retarded Green's function of the 
classical oscillator. The normal ordering is defined primarily for q{t), but, 
since the difference between it and qj{t) is a c-number, it applies equally to 
qj{t). By making use of fl26|) . Eq. fl25|) may be written as 

$R(r/;j) = $ci(r/;j)$in(r?), (27) 

where 

$d(r/; j) = exp J dtr]{t)qj{t) = exp J dtdt'r]{t)D^{t - t')j{t'), (28) 

and 

$in(r/) = ^:exp J dtr^{t)q{ty}j . (29) 

This way, in the response formulation, response properties of the oscillator 
separate from those of the initial state. Moreover, only the latter may be 
quantum. The functional $ci(^?;j) describing the quantum response of the 
oscillator is by itself a fully classical object. It is a characteristic functional for 
products of a classical c-number field emitted by a classical c-number current 
in accordance with laws of classical mechanics: 



qj{ti) ■ ■ ■ qj{tm) 



5r]{ti) ■ --driit^) 



(30) 

r]=0 



In the response formulation, the quantum oscillator allows the most classical 
interpretation possible, assuming that this interpretation remains fully and 
consistently within the laws of quantum mechanics. 



The key result of paper [l[ is that, for the harmonic oscillator, the response 
formulation is formally equivalent with Schwinger's closed-time-loop Green's 
function formulation. In the following fundamental relation between the 
characteristic functionals was derived, 

= ^K{ri,(y + j). (31) 
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where the two pairs of functional variables, ri± and 77, a, are coupled by an 
invertible transformation, 



ri{t) = -i ri+{t) -V-it) 



ait) = h[r]^^\t)+7]'-'{t) 



(32) 
(33) 



Here, g^'^^t) denote the frequency-positive and negative parts of a function 
g(t), defined by dropping the corresponding half (negative or positive) of its 
Fourier-spectrum 



5,(±)(t) = J ^ e-iujte{±uj)g^, = J dt eiujtg{t). 



(34) 



We use the terminology of quantum optics; in the quantum field theory, the 
terms frequency-positive and frequency-negative are often swapped j^. The 
concept of the frequency-positive and negative parts will be used through- 
out this paper; the necessary facts regarding it are summarised in appendix 
lAl Getting back to the characteristic functionals, with j = equation fl3Tl) 
becomes, 



(35) 



For purposes of this paper, the characteristic property of the response formu- 
lation of the quantum oscillator is that it is structurally identical with the 
response picture of a classical oscillator. For the latter, the general solution 
for the displacement is 

q{t) = qUt) + qj{t), (36) 

where qinit) is the free oscillation for j = 0, and qj{t) is defined by fl2^ . If 
the initial condition is random, the classical oscillator is characterised by the 
stochastic moments of the displacement, q(ti) ■ ■ ■ q(tm), with the bar denoting 
the classical stochastic averaging. The characteristic functional for these is 

$R(r7; j) = exp J dtr]{t) [q^{t) + q^it)]. (37) 

Denoting the classical functional fl37|) by $r(?7; j) is perfectly justified, because 
Eq. fl37j) holds equally for the quantum functional fl25l) . For the latter, qinit) 
must be defined as a matrix element over a coherent state 

qin{t) = («|g(t)|«), (38) 
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and the bar should be interpreted as an averaging over the P-function cor- 
responding to the initial state. If the P-function is positive, the classical 
and quantum response pictures are indistinguishable; with a nonpositive P- 
function they remain structurally identical. 



4 The Kubo and Schwinger currents for interacting bosons 



Our first step towards the generalisation of the results for the harmonic oscilla- 
tor to interacting bosons will be finding a replacement for Eq. fl26|) expressing 
the operator "with the source" by the operator "without the source." It goes 
without saying that expecting anything as simple as fl2B]) for the operators of 
interacting fields would be pointless. As was already noted in section 12. in 
principle a connection between Qj(t) and Q(t) is provided by Eqs. (fTT!) and 
ffT2l) . This connection becomes much simpler and more transparent if rewrit- 
ten in terms of the charactestic functionals. For the oscillator, we observe that 
functionals $(?7_,ry+;j) and given by (fT7|l and (ITS]) , respectively, 

are expressed by the same functional $r, hence there must exist a relation 
between them. It is found by noting that the replacement rj ^ ri,a ^ a + j 
amounts to — * ?7± + j/h, so that 

^r]^,T]+;j) = ^(^ri^ + ^-,ri+ + £j. (39) 

Our immediate goal is to prove that this relation holds not only for the har- 
monic oscillator, but also in general for any interacting bosonic system (cf. the 
opening paragraph of section [3]). We note that Eq. (!39|) is sine qua non for 
the very existence of a response formulation of interacting bosons. It makes 
obvious the fact that Q{t) already contains full information on the response 
properties of the system. In a sense, the rest of the paper merely clarifies the 
physical content of Eq. (|39|) . 

We start the proof of equation fl39l) from interpreting $(?/-, ?7+; j) as Schwinger's 
amplitude for the evolution forward and backward in time [6|]: 



$(r/_,r/+;j) = (5l5+), (40) 
where S± are the forward-in-time and backward-in-time S-matrices, 

5±=W±(oo), (41) 
with IA± (t) being forward- in-time and backward- in-time evolution operators. 



U±{t) = T+ exp 



oo 



(42) 
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With 



(43) 



definition (l40l) obviously coincides with (fTTl) . cf. also section 4.3.1 of Ref. pj. 
For simplicity, in Eqs. (l40|) - fH2l) we assumed that the currents j± are real. We 
firstly verify Eq. (139|) for real arguments, then extend this relation to complex 

v±it)=j±it)/h. 

Calling U±{t) evolution operators implies that we can define them as such 
in some dynamical approach. Proof of fl39l) in fact reduces to finding such 
approach and carefully writing down all definitions. Namely, consider the fol- 
lowing pair of Hamiltonians in the Schrodinger picture. 



Tisi (t) = Tt^os {t) +3{t)Qs+ 3± {t)Qs. 



(44) 



With j± = they reduce to the Hamiltonian of the system in Fig. [T] given by 
Eq. ([5]), while j±(t)Qs has been added so as to estabhsh connection with Eqs. 
fH0l) - P3l) . Equation fl39l) will be found as a consequence of the trivial fact, that 
the two constituents of the interaction may be treated either sequentially, or 
concurrently by combining them into a single term 



j{t)Qs+j±{t)Qs= Jit)+J±it) Qs. 



(45) 



Considering the interaction terms sequentially amounts to representing the 
evolution operators related to Hamiltonians flH|) as products of three factors, 

U^{t)U,{t)Mt). (46) 

corresponding to the three terms comprising the Hamiltonians (jUl) . The fac- 
tors lAoit) and Uj{t) were introduced in section [2TT| cf. Eqs. flS]l- ffT^ . while 
the operators U±{t) here are the same as given by P2|) . Indeed, they obey the 
equations 



im±{t)=3±{t)Q,{t)Mt), 



(47) 



solved by fl42l) . Furthermore, considering the interaction terms concurrently 
means replacing lAj{t)U±{t) in fj46l) by Uj±{t), 



U,±{t)=lJ,{t)lJ^{t). 
This pair obey the equations 

iW^j^it) = [j(t) +j±(t)]Q(t)^±(t), 

also solved by T-exponents, 



(4J 



(49) 



Uj±{t) = T+ exp 



dt' 



m+Mt') Qit'] 



(50) 
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Importantly, the T-exponents contain Qj{t) whereas flHIl) contain Q{t). 
By making use of (HHIl then taking t — oo we find 



5± = S] T+ exp ^ -- / dt' j{t') + j±{t') Q{t') \ , 



(51) 



where Sj = Uj{oo) and S± are given by Eq. fllTl) . Equation fl3^ is then 
recovered by substituting (15T|) into (HOj) . using that 



(52) 



and recalhng (Il3 



To extend these considerations to complex j±, it suffices to replace i-{t) — > 
j*_{t) in Eqs. (H^ -(|5T]) while preserving all relations involving j+{t). In par- 
ticular, we then have 



U- (t) = T+ exp 



(53) 



so that 



Wl(t)=T_exp[^ r rft'j_(t')Q.(0 



(54) 



leaving equation fl4UI) unchanged. Importantly, the Kubo current remains real, 
so that none of the physical quantities occuring in Eqs. (Hlll - flHTl) are affected; 
the only affected ones are the evolution operators U±{t) and Uj±{t) which are 
anyway purely formal. Equation (l39l) is thus also proven for complex r]±{t) = 

3±it)/h. 



5 Response formulation of interacting bosons 



5.1 Preliminary considerations 



Response reformulation of interacting bosons is a straigtforward generalisa- 
tion of Eq. fl2^ for the harmonic oscillator. We simply replace the operator 
qj{t) by the Heisenberg field operator Qj{t), and the normal ordering of the 
free-field operators by the time-normal ordering of Heisenberg operators. The 
latter was introduced by Kelly and Kleiner [3] in the context of macroscopic 
photodetection, see also 0,SSIl2|. The "interacting" generalisation of ^ 
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then reads 



°° 1 / /■ \™ / - - \ 

= 1+E^ /^O viti)---vitm){T:Qj{h)---Qj{t^y), (55) 

m=l ^ 



where the notation T: ■ ■ ■ : for the time-normal ordering is borrowed from Man- 
del and Wolf [3]. Within the approximation of slowly varying amplitudes, the 
time-normal averages in ( l55l) may be calculated by expanding each operator 
into the frequency-positive and negative parts, 

Q,{t) = Qf\t)^Qti\t). (56) 
then applying the Glauber-Kelly-Kleiner (GKK) definition. 



T:Qf ^(ti) ■ ■ • Q}-^(t^)Q}-^^(Wi) • • ■ QY\t 



T-Q\-\h) ■ ■ ■ Q\-\t^)T+Q^^\tm+i) ■ ■ ■ Q^^\t^). (57) 



The reason why the GKK definition cannot be used in general will be made 
clear in the next paragraph, and an amended definition will be given in section 
15.31 Discussion of this amended definition will, to a large extent, be subject of 
the rest of the paper. 

Physics behind 0551) is best understood if we consider the alternative way 
of introducing the response formulation: by analogy with the decription of 
the Gedanken experiment in Fig. [1] in the classical stochastic field theory. As 
was mentioned in section [3l the quantum and classical response formulations 
of the harmonic oscillator are structurally identical. They differ only in the 
replacement of the classical stochastic averages by the normal averages. We 
now show that Eq. flS^ follows if we replace the stochastic averages in the 
classical description of the experiment in Fig. [1] by the time-normal averages. 
This means, in particular, that the structural identity between the classical 
and quantum response pictures is a general property of bosonic systems. 

In the most general classical approach, the field seen by the detector in Fig. 
[1] is a c- number stochatic field Qj{t) dependent on the external source jit). 
It is fully characterised by its stochastic moments Qjiti) ■ ■ ■ Qjitm) ■ In turn, 
the field response is fully characterised by the classical stochastic response 
functions 



(5J 

j=0 



As was first shown by Glauber and later extended to Heisenberg fields 
by Kelly and Kleiner [3] (see also 0, S, 0, [l3]), the quantum generalisation of 
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this classical picture consists of replacing the classical stochastic averages by 
the time-normal averages of the quantum field 



Qjih) ■ ■ ■ Qjiu, 



(59) 



The response of the quantum field is then characterised by the quantum- 
statistical response functions 



V 



(m,n) 
R 



1? ■ ■ ■ ; 



f -t' ■ ■ ■ f') 



5^{T:Q^{t,)---Q,{tr, 



5" 



6r]{ti) ■ ■ ■6r]{tm) 



\v=o 



(60) 



j=0 



where $^(77;^) is the characteristic functional of the time-normal averages of 
Qj{t) given by (|55|). 

In addition to clarifying the physical content of the definition fl55|) . Eqs. fl60l) 
show that interpetation of the functional $r(?7; j) is by construction two-fold. 
Viewed as a functional of rj with j a parameter, $R(?7;j) is a characteristic 
functional of the time- normal averages of operator Qj{t). Viewed as a func- 
tional of two variables, $^(17;^) is a characteristic functional of the quantum- 
statistical response functions of the quantum field. 



E 

m+n>l 



mini 



dt 



m+n 



vih)---vit^mt[)---jit'j 



X V 



(m.n) 



it 



1; ■ ■ ■ 5 



tm] ti, - ■ ■ , t„), (61) 



where ( / dt) stands for integration over all ts, cf. also endnote [13|. Con- 
sistency of these two interpretations of <I>r(?7;j) is warranted by equations 



5.2 Failure of the naive generalisation 



A discouraging observation is that the functional $r(?7; j) defined by (l551) - flFn) 
cannot be mapped exactly on $(r7_, j) given by ffTSl) . Consider, e.g., the 
time-normal average of three field operators (T:Qj{ti)Qj{t2)Qj{t3):Y Separa- 
tion of the frequency-positive and negative parts of a function is conveniently 
expressed as an integral operation, 

g(^\t) = Jdt'6^^\t-t')g{t'), (62) 

where 5^^-* (t) are the frequency-positive and negative parts of the (5-function, 
cf. appendix \M On assuming that ti > t2 > t^ and pairwise combing terms 
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we find 



(63) 

The integration here is over all possible values of t[, tg, ^3 including t'2 < t[,t'^. 
However, for this order of arguments, the product Qjit[)Q.j(t2)Qj{t'^) is not 
double-time ordered. Indeed, in a double-time-ordered product, the times are 
arranged into a first-increase-then-decrease sequence, which t[ > < t'^ is ob- 
viously not. One result of paper [ij is that for the free- field operators such 
offending contributions cancel. Unlike the normal products of the free-field 
operators, the time-normal products of the Heisenberg operators as defined 
by Glauber and Kelly and Kleiner cannot be expressed solely by the double- 
time-ordered products of the field operators. 



5.3 Redefining the time-normal ordering 



To find a way of dealing with the complication we have just encountered, let 
us have another look at the results for the harmonic oscillator. Combining 
equations ([2SD and (^11) we find 

:exp j dtrj{t)qj{t): 

= r_exp J dtri^^\t)qj{t)T+exp J dtri^-\t)qj{t). (64) 

Dropping the quantum averaging here is justified, because Eqs. fl25|) and fl3Tl) 
hold for an arbitrary quantum state. The LHS of flM|) is an operator-valued 
characteristic functional of the normal products of qj{t), so that this relation 
may be used as an alternative definition of the normal ordering. We follow 
this pattern and redefine the time-normal ordering of Heisenberg operators as 

TiexpJ dtr]{t)Qj{t): 

= T_ exp J dt r]^+^ (t) Qj (t) T+ exp J dt r]^-\t) Qj (t) . (65) 

By virtue of shorthand ([2]), this defines the time- normal products not just for 
a nonlinear oscillator, but for an arbitrary bosonic field. In the field case, the 
positive and negative-frequency parts always apply to the time variable. 



17 



Unlike (16^ . Eq. (l65i) is not equivalent to the conventional definition by Glauber 
and Kelly and Kleiner (GKK definition). It is therefore instructive to under- 
stand the difference between definitions (1571) and (1651) . There is a very simple 
incorrect way of deriving the GKK definition from Eq. (I^Sj) . by using the 
formula 

J dtr]^^\t)Q,{t) = J dt7]{t)Qf\t). (66) 

This formula follows from Eq. flB^ and from 

5(±)(t-t') = 5W(t'-t), (67) 
cf. Eq. (lA.SP in the appendix. Then, indeed, 



' dt7]^^\t)Qj{t) = J dtdt'rj{t)6^^\t - t')Qj{t') = J dt'7]{t')Qf\t'). (68) 
Applying Eq. (1^ to (1^ we have, 

Tiexp j dtr]{t)Qj{t): 

= T_exp J dtT]{t)Q^'\t)T+exp J dt7]{t)Qf\t). (69) 

It is easy to check that this formula is equivalent to the GKK definition ([5 



To understand the flaw in this derivation is to understand the difference be- 
tween ( 1571) and (IB51) . Consider a term in (IB51) contributing to T:Qj{t)Qj{t'):, 

j dtdt'7]^-\t)r]^-\t')T^Q,{t)Qj{t') 

' dtdt'r]{t)7]{t') [ dt"dt'"6^+\t -t")T+Q,{t")Q,{t'"), (70) 



where use was made of Eq. (1681) . The inner integral in (1701) expresses separation 
of the frequency-positive parts of T^Qj{t")Qj{t"') in respect of t" and t'", cf. 
Eq. (162|) . In Eq. (!70|) (and also in general in (165|) ) the time orderings come 
first, and the separation of the frequency-positive and negative parts second. 
In the GKK definition, separation of the frequency-positive and negative parts 
comes first, and time orderings second. The flaw in the above derivation is that 
it assumes implicitly that these operations commute. In fact they do not, and 
this is exactly how the difference between fl57j) and fl65l) originates. 

It is convenient to preserve the customary notation for the time-normal oper- 
ator products by introducing the shorthand 



T4 



Qf\t)---] = j dt'5'^^\t-t')T^[---Q,{t')---]. (71) 



18 



Using this shorthand means that the GKK definition (157|) has been replaced 

by 

T:Q'j-\h) ■ ■ ■ Q^f\t^)Q^;\t^+i) ■ ■ ■ Q^^\tm): 

= i^j dt'y 6^-\h - t[) ■ ■ ■ 5^-\t^ - t'J5^^\t^^, - ■ ■ ■ - C) 

X TMAt'i) ■ ■ ■ Qi(4)2^+Q.(4+i) • ■ • QjiC). (72) 

The order of the separation of the frequency-positive and negative parts and 
of the time-ordering here is obviously reversed compared to (1571) . 

Usefulness of shorthand (171]) is in fact twofold. Firstly, expressions involving 
this shorthand may be handled algebraically exactly as those with the GKK 
definition. For instance, the formulae 



T:Q,{t): = T:Q'^\t) + Q^-\t): 



T:Q,{t)Q,{t'): = T:[Q'^>it) + Qy>{t)\ [Ql^\t') + Q^f\t')\: (73) 
= T+Q;+)(t)Q(+)(t') + Q5-)(t)Q5+)(t') 

etc., hold with either definition. Secondly, and more importantly, definitions 
(1571) and (172|) need not be distinguished within the quantum-optical paradigm. 
Indeed, analyses in [i| verify that, for the free fields, definitions (157|) and (1721) 
are exactly equivalent, so that there should exist physical conditions under 
which these definitions remain approximately equivalent. It is easy to see that 
one such condition is the approximation of the slowly varying amplitudes. As 
a consequence, our redefinition of the time-normal ordering does not affect 
interpretation of any of the experimental facts in quantum optics involving 
this concept. The reader who is only intererested in applying our results within 
the conventional quantum-optical paradigm may ignore our specifications. The 
reader who is interested in how the definitions initially developed within this 
paradigm must be changed in order to apply to relativistic quantum fields 
(say) should remember that here, strictly speaking, (157|) is just a shorthand 
for (d. 



5.4 Response reformulation of bosons 



The response formulation of interacting bosons is by definition introduced by 
Eqs. ( 1551) . (160|) and (1611) . where the symbol of the time-normal ordering must 
be understood according to ([H5D in place of (lFr|) . The immediate consequence 
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of this definition is that the fundamental relation (l3Ti) connecting the closed- 
time-loop and the response formalisms holds also for the interacting bosons. 
Indeed, on averaging Eq. (^5^ and applying definitions (ITS]) and (jSSD we find 

$R(r/,j) = ^-tr]^+\tr]^-^;j). (74) 

It is easy to see that Eq. (!39l) affords a trivial generalisation, namely, for any 
ji and ja, 

HV-,V+;Ji+J2) = '^(v- + j,V+ + j,J2y (75) 

This relation is verified by applying Eq. fl39|) to both sides of it. Replacing 
J ^ J + cr in Eq. fl7i|) and using Eq. fl75l) then results in 

*r(^,^ + j) = (^-^r/W + ^,ir7^~) + = $ (r^^, 7]+; j) , (76) 

where the pairs of functional variables ?7±(r, t) and //^(r, t), cr^(r, t) are coupled 
by the substitutions 

ri>^ir,t) = -t['ntir,t)-r]'^ir,t)], 

a^(r,t) = ;i[<(+)(r,t)+r7^("^(r,t)], ^ ^ 

^ 1 (78) 
ry^l(r,t) = -ir7^W(r,t) + -fT'^(r,t). 



We have thus recovered relation fl3T|) for interacting fields, with substitutions 
f l32l) . fl33|) replaced by their field versions. Substitutions fl77j) . fl78l) are exactly 
those found in section 3.1 of [l| for noninteracting bosonic fields ensuring 
overall consistency. Note that applying Eq. fl75]) in deriving Eq. fITB]) implies 
that a^{r,t) is real. This assumption is not an obstacle, cf. the discussion in 
section 2 of Ref . jl| . We return to this condition shortly. 

The logic of the above argument may be inverted by using substitutions flTTl) 
and ( 1781) as a starting point in place of Eq. ( I65|) . The functional <I>r(?7; j) is 
then defined a priori by applying (ITSl) to $(r7_,77+). 



$R(r7; a) = <I>(r7_, 77+) , (79) 

while equations (^S^, (IBUl) and (1^ serve as equally a priori definitions of 
the time-normal averages and response functions (independent of the GKK 
definitions). Equation (176|) becomes a consequence of (175!) and (!79|) . while Eq. 
( !7il) is found as a particular case of (176|) with a = 0. The formula for the 
time-normal products (l65l) is then recovered from (17^ by noting that the 
latter must hold for an arbitrary state of the field. Importantly, irrespective 
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of whether we start from Eq. (165|) or from Eq. (1791) . the overall consistency of 
our definitions depends on Eq. (!39|) . 

In essence, we have generalised the results for the oscillator to quantum fields 
by postulating that substitutions fl32|) and fl33l) are always applied to the time 
argument and do not touch field labels. This way of extending the results for 
the oscillator to fields avoids the concept of mode, making our considerations 
applicable, in particular, in the case of arbitrarily strong interactions, when 
introducing modes may be difficult or impossible. 

Despite the generality of our definitions, we will be able to satisfy explicit 
causality in the response formulation and trace its links to Glauber's photode- 
tection theory and Kubo's linear response theory. In particular, we shall show 
that all reality and causality properties characteristic of the classical stochastic 
response functions fl58l) also hold for the quantum-statistical response func- 
tions flBUl) . It is easy to see that the time-normal averages introduced by defini- 
tion (I79p are real quantities; the V^'^^s are therefore also real. This becomes 
obvious if 7], a are chosen real; this choice was in fact already implied when 
deriving Eq. ([76]). Then //_ = r]l making $r(?7; a) = ?7+)|^^=;j-i„±i^(T) 
real. Coefficients of the Taylor expansion of a real functional, cf. Eq. fl55|) . 
are obviously real quantities. Causality in the response formulation will be 
demonstrated in section [71 

It cannot be emphasised strongly enough that equations (!65!l - ([79l) constitute 
a radical deviation from the logic of leading considerations in section 15.11 
Equation fl60p implies that the concept of time-normal ordering is known, and 
uses it to define the response functions, whereas equation (17^ postulates the 
response formulation thus redefining the time-normal ordering. This change 
of the logic, we remind the reader, derives from the fact that the response 
formulation based on the conventional time-normal ordering as defined by 
Glauber and Kelly and Kleiner [3, does not map onto the Green-function 
formulation. This forces us to consider Eq. f[7^ as an a priori definition, 
and the normal averages and the quantum-statistical response functions — as 
purely structural objects defined by expanding functional $r(77; ct) in functional 
Taylor series ([55|1 and (!6T|) . respectively. 

At the same time, response formulation based on definition f[65|) is by construc- 
tion equivalent to the closed-time-loop formalism. The quantum-statistical 
response functions are then expressed by Green's functions f[T3|) . and, impor- 
tantly. Green's functions are expressed by the quantum-statistical response 
functions. In less formal words, everything about the quantum system is in 
its self-radiation and response, and that nothing exists that is decribed by 
quantum mechanics and that cannot be reformulated in self-radiation and/or 
response terms. As has already been mentioned in the introduction, it also 
shows that the information contained in the field operator (and in the intial 
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condition) suffices to describe any scattering experiment performed on the sys- 
tem. It is worthy of stressing that all these results hinge on Eq. (139|) . While by 
itself this formula is in no way limited to the response formulation, this for- 
mulation certainly makes its physical interpretation much more transparent. 



6 Structural response properties of interacting bosonic fields 

6.1 The linear response 

First indication that the way we introduced the response functions makes phys- 
ical sense comes from considering the linear response function V^jl''^\t;t'). It 
is calculated directly from definitions (fT7|) . (!60|) and (179|) . After some straigh- 
forward algebra we obtain 



which is nothing but Kubo's renown formula [3] for the linear response func- 
tion. Details of its derivation based on the general formula for J)^'""^ are given 
in section rr2l Recalling that our definition of the response functions hinges on 
a structural analogy with a linear system, rediscovering Kubo's formula which 
holds for any nonlinear system is quite encouraging. We note also that Kubo's 
formula is the simplest example of the general link between commutators and 
response to be encountered throughout this paper. 

6.2 Quantum response and noncommutativity of operators 

A closer inspection of Eqs. (l39l) and (176|) reveals a fundamental link between 
the operator noncommutativity and the quantum-statictical response of the 
quantum system. Indeed, if we assume that Q{t) for different times commute 
with each other, the time orderings in definitions of ^{ri_,ri^] j) and $(//_, 
may be dropped so we have 







and 
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while (l39ll shows that 



F{r];j)=F{r]). (83) 

Neglecting the noncommutativity of operators thus fully cancels the depen- 
dence on the external source. We further discuss this point in the following 
two paragraphs. 



6.3 General formula for nonlinear quantum- statictical response functions 



The question we address now is how the response functions of the field are 
expressed by the averages of operator Q{t). We note immediately that the very 
fact that such a relation exists is highly nontrivial. Moreover, this relation must 
be inherently quantum. As we demonstrate now, in the quantum field theory 
the response infomation is contained in the commutators. In other words, 
commutators of a Heisenberg field express formally the reaction of the field 
to perturbations. In the classical field theory, the c-number fields commute so 
that this information must be expressed by other means: it enters through the 
Hamilton- Jacobi equations or some other dynamical relations. 



A formula for follows by expanding the RHS of Eq. (1791) into power 

series and comparing them to series fl6T|) . Expansion of the RHS of Eq. fl79l) 
is given by Eq. f|T7|l : by using it we have 



i^i—iV 



X (T^l / dt 



h 



k+l>l 



k\l\ 



-ir]^+\t) 



m tj / dt 



m 

h 



+ ir]^~\t) 




Applying the binomial expansion we have 



k I 



*r(^;j) = i+ E EE 



jk+l(j\l+k 



ktititoto^'- ik-k)\n {i-iy.h'^^^^'^-^ 



X T_ 



dtr]'^+\t)Q{t) 



dtj{t)Q{t) 



X 



dtr]'^~\t)Q{t) 



k-k 



dtj{t)Q{t) 



i-i 



15) 



Comparing this to (!6T]) we see that the terms contributing to D^'""^ are those 
with k + 1 = m, k + I — k — I = n. It is therefore convenient to introduce new 
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summation variables such that 

k = rh, k — k = n, I = m — rh, 1 — 1 = n — n, 

and 

m n ^nl j\n—n 

'^'r(^;j) = 1+ E EE-TT ^Tu ^TT^ 

mT;r>i^o^o "^'("^ -ni)\n\{n- n)\h- 



^6) 



X r„ 



c/tr7W(t)Q(t) 



dt3{t)Q{t) 

dtj{t)Q{t) 



By making use of the shorthand notation ( iTTl) we can rewrite this as 

m n ^n/ j\n—n 

/ /• \ m+n I 

X t/tj r/(ti)---r/(Uj(t'i)---j(0(2^-Q^"nii)---Q^"n^r^) 

X Q(t;)---Q(4)T+QW(W0---Q^+nUQ(4+i)---Q(^^^^ (88) 



where all time integrations are now explicit. By using the definition of the 
response function, cf. Eq. (160|) . we then find 



^(m;n)/, . ,/ 

'^K ) ''■mi ''1) 5 ''n 



n m n 



'-11 



perm m=0 



m!(m — m)!n!(n — n)! 



rij ■ ■ ■ 



X T+QW(r^+0 ■ . • Q(+)(r„)Q(r;^,) ■ ■ ■ QK)). (89) 

Here, ri, ■ • ■ , Xm and t[, - ■ ■ , are, respectively, permutations of ti, ■ ■ ■ ,tm 
and t'l, ■ ■ ■ , t^, and Sperm denotes a summation over a total of m!n! of such per- 
mutations. The average on the RHS of (189!) is symmetric with respect to an ar- 
bitrary permutation within Ti, ■ ■ ■ , r^, within t[, - ■ ■ , r^, within Xm+i, ■ ■ ■ , t^, 
or within r^+i, • • • , r^. Confining the summation to permutations resulting in 
different terms then cancels the factorials in the denominator of the prefactor. 
As a result we have 

.(m;n) , 



l''!' ■ ■ ■ ! '"mi ''1) ■ ■ ■ ; ''n 

i 

h 



n m n 



-t) E E(-if E'(^-2^-nn)---Q(-n^^)Q(rO---Q(r;) 



m=0 n=0 perm ^ 

X T+QW(r^+i) ■ ■ ■ QW(r^)Q«^i) • ■ ■ Q«)\, (90) 
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where I]perm is a summation over all different terms of given structure, i.e., 
such that cannot be transformed one into another merely by permuting factors 
under the time orderings. For n = we find definitions of the time-normal 
averages identical to (175]) . The total number of terms on the RHS of (1901) 
equals 2™"'"" so that expressions for higher order response functions grow very 
bulky. 



Equation fl^U]) does not appear to exhibit any causality properties, and, in fact, 
it takes some effort to prove that such properties hold. Causality conditions 
for pj^'"^ will be discussed in section [TJ where we also derive an alterna- 
tive formula for the response functions making their explicitly causal nature 
evident. 



The overall phase factor in fl90l) depends only on the number of inputs; the 
same applies to the overall power of Planck's constant. The relative signs 
of terms depend on the number of "input" Qs under the T_-ordering. This 
rule can be traced back to the representation of $(//_, as a Schwinger 
amplitude, cf. Eq. (140!) . It is obvious that if we neglect noncommutativity of 
the operators, the sign factor makes all D^™;"-) with n > 1, that is, the response 
functions proper, zeros. This is a direct consequence of equation (!83l) . It clearly 
shows that commutators express response. However it would be incorrect to 
say that the response simply stops functioning in the classical limit h ^ 0. 
Nonzero powers of H are found in the denominators of the quantum formulae 
for the response functions, cf. Eq. fl90l) : in the classical limit, these formulae 

do not give zeros but indeterminates -. Correct classical limit can only be 

achieved by analysing physical details and not by formal means like neglecting 
noncommutativity of operators or setting Planck's constant to zero. 

We illustrate equation fl90l) by constructing a formula for Pj^' . From fl60|) . 



This is the simplest example of a quantum-statitistical response function in 
the true sense of the term. It does not follow from Glauber's photodetection 



> 



j=0 



(91) 
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X 



theory nor from Kubo's linear response theory. For D^' equation (l90l) gives 

X "(T+Q(t;)Q(4)Q(t')) - {mT+mmQ 
\Q{t\)T^m)m) - {t- 

+ c.c. . 

(92) 



In this expression we retain the orderings only where applicable; in the last 
term we put square brackets around Q{t'])Q{t') to delineate the range to which 
the T_-ordering is applied. To emphasise the actual structure of (pOj) . in equa- 
tion ^% we did not use the shorthand (17T]) . 



6.4 Expressing Green's functions by the response functions 



Description of the system in terms of the response functions is by construction 
equivalent to the conventional QFT description in terms of Green's functions. 
By using the inverse substitution (1771) we can write 

$(r/_, r]+) = <l>R(r7, ^)|,=,^__,^^,,=^^(+)+^^(-). (93) 

This way, not only the response functions may be expressed by Green's func- 
tions using (1791) . but also Green's functions may be expressed by the response 
functions using (p3|l . Proceeding in close similarity to how Eq. (pOj) was derived 
we find 

T_Q(ti) ■ ■ ■ Q(4) T+Qit[) ■ ■ ■ Qit[)) = ^ ^(-^)^-^^'-'>+'-^-^ 

fc=o 1=0 

(k+hk+l-k-l)( r- r' ... t^-t(+) ... r'^'K . . r'^'h 

(94) 



X 

perm 



Following the pattern of shorthand (171 p . we introduced here another shorthand 

i,r'(---*»---)^/^«'*'(*-oi'r>(---*--). (95) 

Times and summations in Eq. (!94|) ) have the same meaning as in (pOl) . 

It is instructive to isolate in Eq. (I94p the term with k = k, I = I. This terms 
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equals 

■ ■ ■ ,t,,t[, ■■■X) = {T:Q{ti) ■ ■ ■ Q{tk)Q{t\) ■ ■ ■ my.) . (96) 



This is the only term on the RHS of flMjl that enters with zero power of 
Planck's constant. Other terms are expressed by response functions in the true 
sense of the word which all enter with nonzero powers of Planck's constant. 
This can be written symbolic relation 



T_Q{t,)---Q{t,)T^Q{t[)---Q{t[) 

T:Q{h) ■ ■ ■ Q{tk)Q{t[) ■ ■ ■ Q{t'i):) = hx {Response}. (97) 



The difference on the LHS of Eq. (1971) may be nonzero only because the oper- 
ators do not commute. This equation is another manifestation of the general 
link between noncommutativity and response. 

The terms in the sum denoted symbolically as {Response} may contain addi- 
tional powers of h; some of them may have classical interpretation and others 
may not. Importantly, the hierarchy of powers of h here is not the hierarchy 
of quantum corrections. Indeed, the power of Planck's constant scales simply 
as the the order of nonlinearity. The maximal power of h have the terms with 
k + I — 1 inputs and one output. All these terms are expressed by the same 
nonlinear non- statistical response function 



5^^'-^{Q,{t) 



(9J 

j=0 



where we have used that T:Q{t): = Q{t), cf. Eq. ( 175]) . The function 
may well have a purely classical meaning (for instance, in the mean-field ap- 
proximation). This is just another example of how dangerous it is to formally 
classify quantum effects by the powers of Planck's constant. 

We complete this paragraph with some examples. For k + I = 1 Eq. ([& 
reduces to 



T_Q(t) = new = 2W • (99) 



For A; = 0, / = 2: 

;r+Q(ti)Q(t2))-(r:Q(ti)Q(t2): 

= ih 



V^^-''\t,-A-^)+V^^-^'\t,-4-^)]. (100) 



This is a nonlinear counterpart of Eq. (48) of [1]. This becomes evident if the 
system is homogeneous in time. In this case 

V^^-^'\t,-A-^)=V^^-^'^^'-\t,~t,), (101) 
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and similarly for the second term in fllOOp . For k = I = 1, and assuming time 
homogeneity we find 



ih 



V 



(i;i)(+) 

R 



(ti - h) - V 



R 



(t2-ti) , (102) 



which is a nonlinear counterpart of Eq. (47) of An interesting observations 
is that the RHS and hence the LHS of (11021) remain frequency-positive with 
respect to ti — ^2 also in the general case of an interacting system. 



7 Causality of the quantum response 



7.1 Causality conditions for the response functions 



It is important to realise that nothing in the above warrants causality in 
the response formulation. Formula fl^U]) involves separation of the frequency- 
positive and negative parts of Green's functions flTSl) . In turn this implies 
smearing of the latter from t = —oo to t = +oo. Therefore all terms in fpJl) 
contain a contribution from the future. This "future tail" could well make the 
p(m,n)^ acausal, or only approximately causal. We now show that in fact the 
"future tail" cancels exactly. 

To understand what kind of causality condition is to be expected, we employ 
the structural analogy between the quantum and classical response formula- 
tions discussed in section El The classical response functions fl58l) must vanish 
if at least one input time exceeds all output times. By analogy with the re- 
sponse of a classical stochastic system we formulate the causality condition 
for dI^'") as, 

P|("'"^(ti, ■ ■ ■ ,tm]t\,--- , 4) = 0, max(ti, ■■■ ,tm) < max(t'i, ■ ■ ■ , t'^). 

(103) 

By the same analogy we expect that 

^^r^(;t'i,---,O = 0, n>0. (104) 
This is a quantum analogue of the classical condition 

required by conservation of probability. Indeed, 1 = 1 is a statement that full 
probability is one; (11051) stipulates that full probability does not depend on the 
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external current, in other words, is conserved. In fact, equation fll04p follows 
directly from (!60!) by noting that $r(0; j) = (1) = Trpo = 1 which is full 
probability in quantum mechanics. Equation (11041) thus expresses conservation 
of probability in quantum mechanics. 

We note that a general causality condition for the response functions may 
only be nonrelativistic. Claiming otherwise would mean, for instance, that 
an approach based on the nonrelativistic Schrodinger equation can be made 
relativistically causal merely by rewriting it in response terms. This simple 
example emphasises that relativistic causality is prone to nonrelativistic ap- 
proximations in dynamics and can therefore only be demonstrated within a 
truly relativistic model. 



1.2 Causality of Kubo's linear response 

It is instructive to trace how causality emerges in the linear response function. 
By employing the general formula (I9UI) we have 

Vt\t;t') = -y dt 

X - i) (T+Q{t)Q{t')) - 6^+\t - 1) (Q{t')Q{t)) 

(106) 

If we divide the whole integration region into t > t' and t < t', only the former 
contributes and we find 



+ 6^-\t - t) (Q{i)Q{t')) - 6^-\t - i) {T_Q{t')Q{i) 



- 6^+\t-t) + 6'^-\t-t)\{Q{t')Q{t))\. (107) 



By noting that 



5^^\t-i) +5^-\t-t) = 5{t-t), (108) 



the integration over t may be performed explicitly. Kubo's linear response 
function flHOl) clearly follows. 

This example emphasises the key point: the terms in Eq. (190|) can be pairwise 
combined so that the latest integration time becomes associated with an "in- 
tact" ^-function making the overall order of times meaningful. This idea will 
be used in the general proof of causality given in the next paragraph. 
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1.3 Direct •proof of the causality conditions 



Equations (I103p and (11041) may be traced down to the following simple struc- 
tural property of the double time ordering. Let 11 denote an arbitrary product 
of the field operators where all time arguments are less than some t. Then, 
obviously, 

T_[nQ(t)] = T_[n]Q(t), 
T+[Q(t)n] = Q(t)T+[n]. 

Square brackets here specify the range to which the time orderings are ap- 
plied. The latest field operator in a double-time-ordered product can therefore 
equally be placed under the T+ or T_ ordering: 

T_[n_Q(t)]T+[n+] = T_[n_]T+[Q(t)n+] = T_[n_]Q(t)T+[n+], (no) 

where 11+ and n_ are now two independent operator products with all times 
preceding t. In turn, this means that in the functional Taylor expansion of 
$(?]_, 77+) with respect to ri+,ri_ the latest operator in a double-time-ordered 
operator average always enters in the combination: 

— i 



vAt)-v4t)\ {T4U_]Qit)T4U, 

= v{t)(T-[Tl-m)T4U^]). (Ill) 



so that the latest time is always associated with an "intact" 77 without any 
admixture of a. "Intact" means that the latest rj is not split into the frequency- 
positive and negative parts, making the time order of arguments meaningful. 

To turn these leading considerations into a proof we have to show that 
-ir7+(t) (T_[n_]Q(t)T+[n+]) and zr/_(t) (T_[n_] Q(t) T+[n+]) enter with 
equal weight. We do this by manipulating the range of time integrations in 
(1171) . By making use of the symmetry of the integrand we may confine the 
time integrations to ti < ■ ■ ■ < t^, ^'i < ■ ■ ■ < This cancels the factorials in 
the denominator of the prefactor. We then further split the integration range 
into tm > t'^ and tm < t'^. If tm > t^, we define t = tm, k = m — 1 and / = n. 
The integration region is then 

-00 <ti < ■ ■ ■ <tk <t < 00, 
-00 < t[ < ■ ■ ■ < t'l < t < 00, 

while the prefactor becomes i^~^^{—iy. If tm < t'^, we define t = t'^, k = m 
and I = n—1. The integration region is then again defined by (11121) . while the 
prefactor equals i^{—iy^^. Furthermore, due to the ordering of the integration 
times, the T± orderings may be performed explicitly. The double-time-ordered 
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products for both cases look the same: 

Q(ti)---Q(tfc)Q(t)Q(t;)---Q(t'i; 



;ii3) 



cf. Eq. f lllOp . For every k, I pair, we have two terms with Q{t) coming, respec- 
tively, from the T+ and T_-ordered products. Combining them we get 



k+l+l 



X ■ ■ • Q{U)Q{t)Q{t[) ■ ■ ■ Q{t[)) , (114) 



where the time integrations are confined to (11121) . It is convenient to restore the 
symmetry of time integrations among ti, ■ ■ ■ , tfc and among t[,--- ,ti, resulting 
in 



X T_ 



k+l>0 
t 



kW. 



dt 



dt'r,.{t')Q{t') 



oo 



r7+(t)-r7_(t) 



dt'7]+{t')Q{t') 



>. (115) 



On applying substitution (178|) to this formula we find 
$(r/_,r/+) = l+ J2 m/"^^"^^^^ 



k+l>0 



X ( T_ 



dt' 



h 



xT+l f dt' 



(116) 



Except for one additional time integration and the integration limits, this 
formula is identical to Eq. (18^ and may be manipulated in a similar way, 
resulting in an explicitly causal formula for the response functions, 

^(m+l;n)/, , ,/ ,/ \ 

• \ n m n 

EE(-if^(^-^'i)---^(^-0 



m=0 n=0 



E YT-Q(--)(ri) ■ ■ ■ Q(--)(^^)Q(n) ■ ■ ■ QK) Q{r) 



perm 



X T+Q(+'-)(r^+i) ■ • • Q(+'^)(r^)Q«^i) • • ■ Q«) ), (117) 
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where we have introduced a generahsation of shorthand (17T!) 



— oo 



dt"6'^^\t' ~t")Tj---Q{t")---' . 



(118) 



Other notation in (11171) is as follows: r, ri, • • • , is a permutation of ti, ■ ■ ■ , tm+i, 
r(, • • ■ , is a permutation of t[,- ■ ■ , t^, and the summation is over all differ- 
ent terms of given structure, i.e., such that cannot be transformed one into 
another merely by permuting factors under the time orderings. 

Equation ( \117\\ is a sum of terms such that in each term one output argument is 
bound to be larger than all input arguments, so that causality condition (11031) 
clearly follows. Due to the "time smearing" occuring in (11181) . this selected 
output argument is not bound to be larger than other output arguments, but 
this is irrelevant. We also see that expansion (I116P lacks terms without 77s. 
This proves Eq. (11041) . 

We conclude this paragraph by constructing an explicitly causal formula for 
V^'^^ alternative to fl92l) . For D^'^^ equation (11171) gives 



where square brackets delineate the range to which the orderings are applied. 
To give more emphasis to the actual structure of (11171) we did not use here 
the shorthand (1118^ . 



8 Conclusion 

We have introduced a response formulation of an arbitrary interacting bosonic 
field in terms of dependence of the time-normal averages of the field opera- 
tor on the external source. While being, by construction, equivalent to the 
conventional Green-function techniques of the quantum field theory, this for- 
mulation exhibits key properties of a classical response picture such as reality 
and explicit causality. Validity of these results outside the quantum-optical 
paradigm requires an amendment to the conventinal definition of the time- 
normal operator ordering by Glauber and Kelly and Kleiner. 
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A The positive and negative frequencies 

The operation of separating the frequency-positive and negative parts of a 
function was defined by Eq. (IMIl . Here we discuss some properties of this 
operation. Obviously, 

g^^\t)+g^-\t)=g{t). (Al) 
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This relation, however, imphes that (7^ is smooth at = 0. An important 
property of is the connection between the frequency-positive and negative 
parts and the convolution: 

J dt'h^^\t - t')g{t') = J dt'h{t - t')g'^^\t') 

" dt'h^^\t-t')g^^\t'). (A.2) 



It is proved by noting that the Fourier-transform of any of the three expressions 
is 9{±u))g^hu}. With h{t) = 6(t) this conveniently expresses separation of the 
frequency-positive and negative parts of a function as an integral operation, 

^(±)(t) = J dt' 6{t-t')g^^\t') = J dt'6'-^\t-t')g{t'), (A.3) 
where 5^'^\t) are the frequency-positive and negative parts of the 5-function, 

S^^\t)= J; (A.4) 

Obviously, 

5(±)(t) =5W(_t), (A.5) 

so that 

' dt6^^\t-t')h{t) = I dt6^^\t' -t)h{t) = h^^\t'). (A.6) 



This formula is instrumental in deriving explicit expressions for the response 
functions; it also shows that 

/ dth{t)g^^\t) = f dtdt' h{t)5'^^\t - t')g{t') = f dth^^\t)g{t). (A.7) 



Notation flM|) becomes ambiguous for complicated expressions and composite 
functions like g{h{t)). Taking the frequency-positive and negative parts of an 
arbitrary expression with respect to t will therefore be denoted as [■ ■ 'J^ ^• 
This way, 

[gimt^ = I dt'5^^\t - t')g{h{t')), (A.8) 

as opposed to 

g^^\h{t)) = J dt'5^^\h{t) - t')g{t'). (A.9) 
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In general, taking the frequency-positive and negative parts of something im- 
phes [■ ■ ■jP''. This is consistent with because, obviously, 

bWf^=^?^^Ht). (A.io) 

Both time inversion and complex conjugation replace e —iut —>■ e iut; hence 
g^^\—t) and g^^^*{t) are frequency-negative while g^^\—t) and g^^^*{t) are 
frequency-positive. By making use of flMl) it is easy to prove that 

[9i't)i^^ = g'^^H-t), [9*{t)i^^ = 9'^^^*{t). (A.ll) 
These relations hold for arbitrary g{t). For an even g{t) = g{—t) we then have, 

^W(-t)=^W(t), (A.12) 

while for a real g(t) = g*(t), 

gi±h^t)=g^^\t). (A.13) 
B Response of charged bosons 

Extention of the results to charged bosonic fields reduces in essence to dou- 
bling all arguments in the relations derived for neutral bosons by applying 
the shorthands fl22|) and fl23l) . The "charged" analog of $(?7_,?7+) is given by 
flMl) . A similar formula with Tj,!pj in place of T^T"^ defines the analog of 
$(r7_, j), while the "charged" analog of fl5^ reads 

r7_,r/+,r/+;j,j*) = $ (^r/_ + |,r7_ + + + (B.l) 

Definition of $r(?7; o") turns into 

<l>R(r/, r^; a, cr) = <l>(77_, r/_, 77+, r7+)|..., (B.2) 
where the ellipsis stands for substitutions (1751) supplemented by 

<(r,t) =2r/'^(-)(r,t) + V(r,t), 

% (B.3) 

ff_{r,t) = -«r/^W(r,t) + -a'^(r,t), 

n 

which is exactly the set of substitutions found in section 3.2 of [H for the 
charged noninterating bosonic field. Equation fl76l) becomes 

^{v-,v-,v+,v+;3J*)\- = ^R{v,v;(^ + j,^ + f)- (b.4) 
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Definitions of tlie time-normal averages and response functions of charged 
bosons take the form 



V 



{m,n;k,l) 
R 



(^1; ' ' ' 1 tmi tl, ■ ■ ■ , tn, t-^, " " " , tj^, t^, ■ ■ ■ , 



T :P^{ti) ■ ■ ■ P,{tm)P]{h) ■ ■ ■ P]{tn): 



(B.5) 



driiti) ■ ■ ■ 6T]{tm)Sr]{ti) ■ ■ ■ 6ri{tr 



v=v=i=j*=o 



cf. Eq. f l60l) . The formula for 



{m,n\k,l) 



derived in full analogy to fl^ . looks 



V 



(m,n;k,l) 
R 



• \ k+l k I m 11 _ _ I 

-|) E E E E(-i)'^T ■ --^F^-Kr. 

A:=0 i=0 ™=0 «=0 perm \ 



X P^-\fn+,) ■ ■■P^-\fn)P{f{^,) ■ ■■Pin)). (B.6) 

Here, ri, ■ ■ ■ , r^, ri, ■ ■ • , r„, t{, ■ • • , and • • ■ , r/ are permutations of, re- 
spectively, ti, ■ ■ ■ , tm, ^1, ■ ■ ■ ,in, t[,--- ,t'f., - ■ ■ , i'l; ^'perm ^ summatiou 
over permutations resulting in different terms on the RHS of flB.6l) . Expres- 
sions for the time-normal averages of charged fields follow from flB.6l) as a 
special case k = I = 0. The response functions of charged bosons obey the 
causality properties. 



/j-v(m,n;fc,/) / , 
^R l^^l' 



f -f 



' ''fci '-1' 



max{ts,ts} < max{t'g,t'g}, 



(B.7) 



vS^'-^''^\;;t[,...,t',;i[,...,i\) = 0, (B.8) 
and the reality property, 

^R \tlj ' ' ' 5 tm] ii, ■ ■ ■ , in] ti, ■ ■ ■ , tj^] ii, ■ ■ ■ , ii) 

= T)^ ' ' ' \ti,--- , t„; ti, ■ ■ ■ , ■ ■ ■ 5 ■ ■ ■ }tk)- (B-9) 

Similar to how reality and causality properties of neutral bosons are those of a 
real c-number field, properties flB.7p - flB.9l) of charged bosons are clearly those 
of a complex stochastic c-number field. 
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